Search for: All records

Abstract We perform the numerical simulations of axisymmetric, relativistic, optically thin jets under the influence of the radiation field of an accretion disk. We show that starting from a very low injection velocity at the base, jets can be accelerated to relativistic terminal speeds when traveling through the radiation field. The jet gains momentum through the interaction with the radiation field. We use a relativistic equation of state for multispecies plasma, which self-consistently calculates the adiabatic index for the jet material. All the jet solutions obtained are transonic in nature. In addition to the acceleration of the jet to relativistic speeds, our results show that the radiation field also acts as a collimating agent. The jets remain well collimated under the effect of radiation pressure. We also show that if the jet starts with a rotational velocity, the radiation field will reduce the angular momentum of the jet beam. 

GW170817 and GRB 170817A provided direct evidence that binary neutron star (NSNS) mergers can produce short gamma-ray bursts (sGRBs). However, questions remain about the nature of the central engine. Depending on the mass, the remnant from a NSNS merger may promptly collapse to a black hole (BH), form a hypermassive neutron star (HMNS) which undergoes a delayed collapse to a BH, a supramassive neutron star (SMNS) with a much longer lifetime, or an indefinitely stable NS with a mass below the TOV limit. There is strong evidence that a BH with an accretion disk can launch a sGRB-compatible jet via the Blandford-Znajek mechanism, but whether a supramassive star can do the same is less clear. We have performed general relativistic magnetohydrodynamics simulations of the merger of both irrotational and spinning, equal-mass NSNSs constructed from a piecewise polytropic representation of the nuclear SLy equation of state, with a range of gravitational masses that yield remnants with mass above and below the supramassive limit. Each NS is endowed with a dipolar magnetic field extending from the interior into the exterior, as in a radio pulsar. We examine cases with different initial binary masses, including a case which produces a HMNS which collapses to a BH, and lower mass binaries that produce SMNS remnants. We find similar jetlike structures (helical magnetic field structures, a magnetically dominated evacuated funnel, and mildly relativistic outflow from the poles) for both the SMNS and HMNS remnants that meet our basic criteria for an incipient jet. The outflow for the HMNS case is consistent with a Blandford-Znajek (BZ) jet. There is sufficient evidence that such BZ-powered outflows can break out and produce ultrarelativistic jets so that we can describe the HMNS system as a sGRB progenitor. However, the incipient jets from the SMNS remnants have much more baryon pollution and we see indications of inefficient outflow acceleration and mixing with the surrounding debris torus. Therefore, we cannot conclude that outflows from SMNSs are the progenitors of sGRBs. 

Every numerical general relativistic investigation starts from the solution of the initial value equations at a given time. Astrophysically relevant initial values for different systems lead to distinct sets of equations that obey specific assumptions tied to the particular problem. Therefore, a robust and efficient solver for a variety of strongly gravitating sources is needed. In this work, we present the OpenMP version of the Compact Object CALculator (COCAL) on shared memory processors. We performed extensive profiling of the core COCAL modules in order to identify bottlenecks in efficiency, which we addressed. Using modest resources, the new parallel code achieves speedups of approximately one order of magnitude relative to the original serial COCAL code, which is crucial for parameter studies of computationally expensive systems such as magnetized neutron stars, as well as its further development towards more realistic scenarios. As a novel example of our new code, we compute a binary quark system where each companion has a dimensionless spin of 0.43 aligned with the orbital angular momentum. 

Triaxial neutron stars can be sources of continuous gravitational radiation detectable by ground-based interferometers. The amplitude of the emitted gravitational wave can be greatly affected by the state of the hydrodynamical fluid flow inside the neutron star. In this work, we examine the most triaxial models along two sequences of constant rest mass, confirming their dynamical stability. We also study the response of a triaxial figure of quasiequilibrium under a variety of perturbations that lead to different fluid flows. Starting from the general relativistic compressible analog of the Newtonian Jacobi ellipsoid, we perform simulations of Dedekind-type flows. We find that in some cases the triaxial neutron star resembles a Riemann-S-type ellipsoid with minor rotation and gravitational wave emission as it evolves towards axisymmetry. The present results highlight the importance of understanding the fluid flow in the interior of a neutron star in terms of its gravitational wave content.